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Abstract 

We present a new one parameter family of second derivative discontinuous solutions to the 
simplest scale invariant linear ordinary differential equation. We also point out how the construction 
could be extended to generate families of higher derivative discontinuous solutions as well. The 
discontinuity can occur only for a subset of even order derivatives, viz. ,2nd, 4th, 8th, 16th, ....The 
solutions are shown to break the discrete parity (reflection) symmetry of the underlying equation. 
These results are expected to gain significance in the contemporary search of a new dynamical 
principle for understanding complex phenomena in Nature. 
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1 Introduction 



Let us consider the simplest scale invariant linear ordinary differential equation (ODE) 

dr ,„ > 

A well known fact of the standard analysis (calculus) is that this equation , and more generally, an 
ODE of the form % = f(r, t) possesses a unique C r solution in a neighbourhood of (to, To), tq = r(to) 
when / £ C r . This is the Picard's theorem. The standard solution of eq([T]), as admitted by the 
Picard's theorem, is of course T s (t) = t when the initial condition is chosen as r(l) = 1. Clearly, this is 
a C°° solution. In the following section, we, however, present, in contradiction to the Picard's theorem, 
a new exact solution of eqflT]), which is only C 1 . The new solution breaks the reflection symmetry 
(t — > —t) of the underlying ODE spontaneously. We interprete the solution in Sec. 3 and show how 
this defines a nonstandard extension of the real number system. Consequently, every real number is 
identified with an equivalence class of a continuum of new, infinitesimally separated elements, which 
are in a state of random fluctuations. We also point out how the new construction generalizes to 
create higher derivative (4th, 8th,...) discontinuous families of solutions to the above equation. The 
results presented here are new and sums up our earlier investigations in a mathematically 

rigorous manner. These rigorous results should become important in the contemporary search of a new 
dynamical principle for understanding complexity at different levels of life viz., physical, biological, 
financial and social [1, 0, IE)] , besides their deep significance in mathematical analysis. RefQ] contains 
some preliminary results and definitions indicating their role in defining a new framework for addressing 
living systems (see also @]). 

2 C 1 solution 

Let us first demonstrate the existence of such a solution in the neighbourhood of t = 1. The solution 
would only have a continuous first order derivative at t = 1. We begin by defining our notations. 

Let t n ± = 1 ± %, t Q = t, < ry n << 1, a n = 1 + e n , n = 0, 1, 2, . . . , and e = 0, < e n < l(n ^ 0), 
such that e n — > 0, asn — > oo (we retain the symbol «o f° r the sake of symmetry). Next, we write 
t' n ± = 1 ± a n rj' n , so that a n t n _ = t' n _. Consequently, rj' n = r\ n — Here, a n (and e n ) are scaling 
parameters. A useful choice, however, is e n = e 2 ™, e = e\ (the reason will become clear later). As will 
become evident, r) n+ \ = a^'n- 

To construct a nontrivial solution (with the initial condition r(l) = 1), we begin with an initial 
approximate solution in the small scale variable 770, and then obtain recursively self-similar corrections 
over smaller and smaller scale variables t]q, t/q, To this end, let 

rit) = { l~ !| 1 1 J , r_(t_) = (l/t' + )f^( Vo ), T+ (t+) = U (2) 

be an exact solution of eq(pQ). This is obviously true for the right hand component t + . To verify 
the same for the nontrivial component t_, we differentiate it with respect to i_, and use the scale 
invariance of eq([I]). Utilizing ctot- = t'_, one obtains 

where /{_ = -jjjijS Vo = a o??o- Consequently, eq© would be an exact solution if and only if /]_ solves 
exactly the self-similar equation 

^-A- (4) 

in the smaller logarithmic variable rn£^"_, where t\- = 1 — a 2 ,^ 2 ^ ^ ~ Vi- The self-similar replica 
eq(@J follows from the equality 



so that T- is an exact solution of eq(pQ). The exact (nontrivial part of the ) solution could thus be 
written recursively as 



r- = C^--j-...-^—f n Mn) (6) 

l + l l+ (n-l)+ 

where f n satisfies the nth generation self-similar equation 

^n— "T7 = fn— (?) 
Oln— 

and t n - = 1 — a^_ 1 r/^ 2 _ 1 = 1 — rj n . We also note that t' + = t + , since ao = 1. 

Plugging in the initial condition t± = 1 at t± = 1 (viz., % = 0), one obtains finally the desired 
solution as 

T _ = C±-}--}-..., r + = t + (8) 

r + r l+ ^2+ 

where C = (0)^4.(0) • • •• Notice that 1, since rj[ = — ei/ai, r/ 2 = ef — 62/0:2, etc, when 770 = 0. 
A remark is in order here. 

The solution eq([8j) follows from eq© only if the sequence {As-} is convergent. In fact, we prove 
that /oo = linin^oo / n _ = 1. Let r n = j--rr— ■ ■ ■ tt— • Then for 770 sufficiently small and e n — > for 

n — > 00, the sequence {r n } is convergent (to a nonzero value), since i' n+ — )• 1 as n — > 00. Accordingly, 
for e > 0, 3iVi such that |r m — r n \ < e for m, n > JVi (m > n). As a result, < h\ < r n < &2) ^2 ~ 
0(1), for n > N2 for a sufficiently large N%. Again, / n _, being defined by eq®, is uniformly bounded 
in a neighbourhood of t = 1, so that |/ n -| < & for n > N2- The desired convergence now follows from 
the Cauchy convergence criterion, since |/ n _ — f m -\ = \t~ 1 \\f m -\\T m — T n \ < k^ l ke Vm,n > N, N= 
max(iV"i, N2). Finally, eq©, in the asymptotic limit n — > 00, yields /oo = A J t n ~ [oo = ^j|t=i = t(1) = 1. 

Now to test the continuity of the derivatives of the solution ([8]) at t± = 1, i.e., at r/o = 0, we note 
that rj n is a polynomial in r/o, of degree 2 n , being defined recursively by rj n = r\ n — ^- , r\ n = c? n _ x rj^_ x . 

As a result ^j- = 0, but 7^ 0, at 770 = 0. One thus obtains 

dr_ 1 gi dr/^ a 2 dr/ 2 

1 — = — m i V \~V~a t)~a ^ \~T~i T)~a h '--J ^1 

dr/ 1 + Vo 1 + ocirfx dr/ 1 + a 2 r/ 2 dr/ 

so that = 1 = Tjp 1 at r/o = which means that the first derivative of the solution is indeed 
continuous for all r/o- However, as is verified easily from eq([9|), the second derivative of r_ at r/o = 
is not zero, as one expects on the basis of the standard solution r s = t. Indeed, one can verify that 
= 2(1 — j^ — . . .) ^ at r/o = 0, unless e n = 0, for all n. In this special case, i.e., when e n = 0, Vn, 

our solution ([8|) reduces to the standard solution, since t_ = ^j- 1 a - 1 4 . . . = 1 — 7/0 = t_. 

It thus follows that the solution ([8]), with nonzero scaling parameters, is indeed nontrivial, because 
of this second derivative discontinuity at r/o = 0, that is at t = 1. In fact, the scaling invariance 
of eq([T|) tells also that, t = 1 could be realized as t — > t/io = 1) so that the nontrivial solution (|8|) 
actually holds in the neighbourhood of every real number to, the 2nd derivative being discontinuous 
at t = to- This is our main result of this paper. Let us note here that r_ = r s _ (1 + 0(770)), besides the 
arbitrariness of the scaling parameters e n . Combining the standard and the new solutions together, 
one can finally write down a more general one parameter class of solutions 

T B (t) = t(l + <f>(t)), 4>(t) = et-\{t) (10) 

Note that 



dt 

because r is an exact solution of eq(pQ) (we fix e n = e 2 ™, on anticipation of Sec. 3). The 2nd derivative 
discontinuity of r, however, tells that can not be considered simply as an ordinary constant [H, 0). 

2.1 Symmetry breaking 

The solution ([8]) is not parity (reflection) symmetric. Let P : Pt± = t T denote the reflection trans- 
formation near t = 1 (Pij = —n near 77 = 0). Clearly, eq([T]) is parity symmetric. So is the standard 
solution t s ± = t± (since Pt s = r s ). However, the solution ([8]) breaks this discrete symmetry spon- 
taneously: rf = Pt + = i_, t± = Pt- = C-f--J—J— . . ., which is of course a solution of eq©, but 

clearly differs from the original solution, t± ^ t±. The explicit breaking of the parity symmetry is 
an important new feature of the solution ([8]). Its significance in dynamical systems would be studied 
elsewhere. 



3 Interpretation 

A basic assumption in the framework of the standard calculus is that a real variable t changes by 
linear translation only. Further, t assumes (attains) every real number exactly. However, in every 
computational problem within a well specified error bar, a real number is determined only upto a 
finite degree of accuracy eo ,say. Suppose, for example, in a computation, 1 is determined upto an 
accuracy of ±0.01, so that 1 here effectively stands for the set l e = {1 ± e, e < eo = 0.01} with 
cardinality c, the cardinality of the continuum. We call e an 'infinitesimally small' real number. 
Now, any laboratory computational problem (experiment) is run only over a finite time span, and 
the influence of such infinitesimally small e's, being insignificantly small, could in fact be disregarded. 
However, if the experiment is allowed to run over an 'infinitely' large period of time, analogous to 
a natural (or biological/ financial) process, such as the fluctuations of a river (ocean) water level 
(heart beats of a mammal/ stock market prices ), which is known to occur over many time scales, the 
contributions from the infinitesimally small scales e need not remain negligible. 

Our solution (|SJ) captures the essence of the above scenario in an elegant manner. Notice that 
our initial ansatz t_ ~ approximates the standard solution r s _ = t- upto 0(r/ 2 ) (i.e., when 

r)o ^ 0, but (O(?7o) = 0). At this level of approximation l e = 1 for e ;$ rf. However, at the second 
level of iteration with the variable ii_ = 1 — 775% the effect of e, being of the same order of t? 2 , can 
not be negligible. Recalling that 1 here should be identified as lie, an appropriate variable can 
therefore be written as at\- = 1 ± e — arffi, which in turn is written equivalently as the rescaled 
variable t'i_ = 1 — where rj^ = rj^^fe/a so that = aii_, a = lie. This scaling transformation, 
which is used recursively in ([8]), thus not only reveals an intrinsically approximate nature of l e , but 
also provides an window to probe the fine scale structure of l e by recursively approaching more and 
more accurately the 'exact' value 1 through the finer and finer scales e 2 , e 4 , . . . successively. ( This 
explains the choice e n = e 2 ".) The exact value of 1 (and, hence, of any real number), as endowed 
in the framework of standard classical analysis, is, however, illusory, because of the fact that the 
solution (JHJ), being 2nd derivative discontinuous at t = 1, is a new exact solution of eq([I]). Further, 
the value of a small real number 770 is uncertain not only upto 0(r/ 2 ), but also because of the arbitrary 
parameter e. We note that two solutions r g and r s are indistinguishable for t ~ 0(1) and t?q << 1. 
However, for a sufficiently large t ~ 0(e^ 1 )(= 0(r?Q~ 2 )), the behaviours of two solutions would clearly 
be different. Moreover, because of the intrinsic (irreducible) uncertainty the new solution r g could 
effectively be written as r g : r_ ~ l/i+, t + = t+, without any loss of generality (and information) at 
the order r/Q ^ 0, t]q = 0. The correction factors in eq© only reveals the approximate (statistical) self 
similarity of the solution over higher order scales r/ 2 , 77 4 , . . . and so on. 

To summarise, the nature of the solution ([8]) depends crucially on the nontrivial scaling a n £_ = t'_, 
as defined above, which shifts 7/0 = at the level of first iteration to 77^ = for the second iteration 



and so on recursively to finer and finer scales. The origin of the second derivative discontinuity is 
obviously an effect of this receding of the zero (0) down finer scales. Let us also point out that the 
scaling crf_ = t'_, does not mean at + = t' + (in fact, t' + = (1 — e)(l + jz^Vo) a *+)> which in turn 
leads to the nontrivial feature that is revealed in the solution (jHJ). 

Now to state the above more formally, we note that the standard solution t s = t defines the 1-1 
identity map R — > R of the real number set R. The generalized solution (|10p now provides a one 
parameter extension of the identity map l e : R — > R. Accordingly, every real number t is mapped to 
the corresponding 'fat' real number t e = t(l + (ft), which is the set of 'infinitely many' new elements 
(with cardinality 2 C ) separated by infinitesimally small numbers (distances) [l|, since (ft is infinitesimally 
small. Accordingly, the set of real numbers R is considered to be an enlarged nonstandard set 0] R 
with infinitesimally small elements, which remain dormant in the standard analysis (i.e., in the context 
of the standard solution t s of eq([T]) ) . This extension of R to the nonstandard set R via the generalized 
solution of the linear ODE eq([T]), is analogous to the extension of the ring of integers to the field of 
rational numbers as the solution space of the linear equation ax + b = 0, a, b being integers. Let us 
note further that eqQ written as 

1 1 + % 1 + air/i 1 + a 2 v' 2 

tells that the variables i+, t' 1+ , . . . behave as independent variables. Writing tC = IIq ^^, a fat real 
variable in the neighbourhood of l e , eq (!12p restates the fact that the solution eq([8]) is indeed an 'exact' 
solution of the equation dl *^ = r in R. In R, however, we have an irreducible 'unknown' component 
f n - ~0(1), which is exactly determined only in the limit n — > oo. Ironically, however, even in this 
asymptotic limit, one fails to improve (determine) the solution with an unlimited accuracy, because 
of the indeterminate parameter e, as reflected in the 2nd derivative discontinuity. The fat variable 
reduces to the ordinary real variable t_ only in the approximation when t]q = (r]o ^ 0). 

The solution ([8]) also tells that a fat real variable t = tl in R can change not only by linear 
translations but also by inversions: i_ — > tZ 1 = i+, in the neighbourhood of t = 1. Clearly, the 
inversion is realized 'exactly' in R, when an (first order) infinitesimal (ft is defined by (ft ^ 0, but 
O(0 2 )=O. Consequently, the change (increment) of a real variable could be visualized as an SL(2,R) 
group action, viz., a combination of linear translation ( for ordinary real variable t following the 
standard solution of eq([I|) and inversion ( for t± G 1). As note already, because of the arbitrariness in 
e, as well as being defined upto 0((ft 2 = 0), the infinitesimal elements (ft enjoy an element of randomness, 
which in turn renders t random. The ordinary non-random real variable t is thus retrieved only under 
an approximation, that is, when (ft = 0. 

Let us remark finally that one can generate 4th or higher derivative discontinuous solutions of 
eq([TJ by introducing infinitesimal scalings at an appropriate level of iteration. Note that the 2nd 
derivative discontinuous solution is obtained when the nontrivial scaling viz., t\- — > t'^_ is introduced 
at the 1st iteration (the ansatz in eq© stands for the zeroth iteration ). Thus, the 4th derivative 
discontinuous solution is obtained when the rescaling of the appropriate variable is used at the 2nd 
level of the iteration, viz., t 2 - — > t' 2 __, giving rise to yet another class of solutions 

T « = C^-^-^-%ldots-^—f n ^r,' n ) (13) 

C+ tl+ l 2 + \n-l) + 

where h+ = 1 + ryg, t' 2+ = 1 - e + (1 + e)^, t' 3+ = 1 - e 2 + (1 + e 2 )[(l + e)r^ - e] 2 , . . ., and so on 
for higher derivative discontinuous solutions. Note that T m = r s (l+0(^)) for an + e < 0(t#). 
Consequently, ln(rl n ^/r s ) ~ t]q" corresponds to the nth order infinitesimals in R. 



4 Conclusion 



We have presented new families of higher derivative discontinuous solutions of the ODE ([T]), which 
do not respect the Picard's theorem. The origin of this violation could be traced to the fact that a 



variable in R may undergo changes (increments) via the extended SL(2, R)-like group actions. These 
solutions break explicitly the parity symmetry of the underlying ODE. The derivation of such solutions 
from any linear first order ODE is obvious. Higher order equations will be considered elsewhere. We 
close with the observation that the results presented here may as well be considered to offer, so to 
speak, an exact proof of the fact that mathematics is inexact. 
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